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Abstract: 

We present the calculation of the most important electroweak corrections to dijet 
production at the LHC and the Tevatron, comprising tree-level effects of O (a s a, a 2 ) and 
weak loop corrections of O (a 2 a). Although negligible for integrated cross sections, these 
corrections can reach 10—20% in the TeV range for transverse jet momenta ky. Our 
detailed discussion of numerical results comprises distributions in the dijet invariant mass 
and in the transverse momenta of the leading and subleading jets. We find that the weak 
loop corrections amount to about —12% and —10% for leading jets with A; T ~ 3 TeV 
at the 14 TeV LHC and &>r ~ 800 GeV at the Tevatron, respectively. The electroweak 
tree-level contributions are of the same generic size and typically positive at the LHC 
and negative at the Tevatron at high energy scales. Generally the corrections to the dijet 
invariant mass distributions are smaller by at least a factor of two as compared to the 
corresponding reach in the &t distributions, because unlike the fcx spectra the invariant- 
mass distributions are not dominated by the Sudakov regime at high energy scales. 
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1 Introduction 



The unprecedented energy regime that is accessible at the LHC allows for the investig- 
ation of the laws of physics at the smallest distances. The inclusive production of two jets 
(dijets) at the LHC, pp — > jj + X, allows for a detailed study of QCD at TeV energies. 
Furthermore, several extensions of the Standard Model predict new heavy particles which 
might be visible via dijet signatures in the detector [lj. Some examples are excited states 
of composite quarks q*, string resonances, new heavy gauge bosons W, Z', etc. Inclusive 
jet and dijet production has been analyzed by the ATLAS 0] and CMS [i[ collaborations 
at a centre-of-mass (CM) energy of 7 TeV giving sensitivity to dijet invariant masses of 
up to 5 TeV andjet transverse momenta of up to 2 TeV at the LHC. The Tevatron ex- 
periments CDF Q and DO j^j have investigated jet production up to transverse momenta 
of several hundreds of GeV. At the current level of experimental and theoretical accuracy, 
the SM is able to describe data quite well. At the LHC design CM energy of 14 TeV, 
the energy reach will even go deeper into the TeV range, so that theoretical predictions 
especially have to carefully include radiative corrections that are sensitive to high scales. 

The results for the production of two jets at leading-order (LO) accuracy in QCD 
had been available jfj long before higher-order corrections were established. Later, in the 
1990s, the differential cross sections to inclusive single-jet and two-jet production were 
discussed at next-to- leading order (NLO) accuracy in QCD 0^9[- Currently, enormous 
effort is put into the calculation of the NNLO QCD corrections to dijet production (see e.g. 



Refs. 10Nl6| and references therein). The purely weak corrections of O (a 2 a) have been 



calculated for the single-jet-inclusive cross section in Ref. [17] . and pre liminary results of 



the weak corrections to the dijet production were published in Ref. [18]. The two results 



however, do not seem compatible with each other. Electroweak corrections also were 



calculated for the related process of bottom-jet production [19 



In spite of their suppression by the small value of the coupling constant a, the elec- 
troweak (EW) corrections can become large in the high-energy domain 0-0. This is 
due to the appearance of Sudakov-type and other high-energy logarithms that result from 
the virtual exchange of soft or collinear massive weak gauge bosons. The leading term is 
given by ct w ln 2 (Q 2 /M{y), where Q denotes a typical energy scale of the hard-scattering 
reaction, M w is the W-boson mass, and a w = a/s 2 ^ = e 2 /{i.ns 2 v ) is derived from the 
577(2) gauge coupling e/s w with s w denoting the sine of the weak mixing angle 6 W . In the 
case of massless gauge bosons, e.g. in QED or QCD, these logarithms are connected to the 
well-known infrared divergences and are cancelled against the corresponding real-emission 
corrections. For the massive gauge bosons W and Z, no such singularities occur, since 
their masses provide a physical cut-off and the additional radiation of real W or Z bo- 
sons can be experimentally reconstructed to a large extent, so that W/Z bremsstrahlung 
corresponds to a different class of processes. Thus, at high scales \Q 2 \ ^> M^, which 
are accessible at the LHC and the Tevatron, the above Sudakov-type logarithms can pro- 
duce large negative corrections, as only some fractions are compensated by unresolved 
W/Z emission 28] . It turns out that large compensations can occur between different 



electroweak logarithms [29j, [30], so that a full fixed-order calculation desirable. 

As it will be discussed in more detail below, a gauge-invariant classification of the EW 
corrections into photonic and purely weak corrections is feasible for our process. Guided 
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by the logarithmic enhancements, we restrict ourselves to the calculation of the purely 
weak corrections in this paper. The calculation can be complemented by the photonic 
corrections at a later time to produce results for the full EW corrections at the order a^a. 

The paper is organized as follows: In Section [2] we set up our conventions (Sect. 12. ip 
and discuss our strategy for the calculation of the NLO corrections (Sect. I2.2j) . The 
numerical results are presented in Section |3l which comprises integrated cross sections as 
well as differential distributions and the comparison to other work. Section H] contains 
our conclusions. 

2 Dijet production in hadronic collisions 

2.1 Conventions and calculational setup 

We consider the hadronic process 

A(p A ) + B(p B ) j(h) + j(k 2 ) + X, (2.1) 

where the assignment of the four-momenta to the respective particles is indicated in 
parentheses. We further assume that the momenta k\ and k 2 are sorted in a descending 
order with respect to their transverse momenta, i.e. &t,i > &t,2> referring to the associated 
jets as the leading and subleading jet, respectively. The hadronic cross section is given by 
the incoherent sum over the different partonic subprocesses that contribute to the final 
state under consideration, convoluted with the respective parton distribution functions 
(PDFs), 

0~Ab(PA,Pb) = ^2 dx b fa\A(x a ,f4)fb\B(Xb,f4)^ab(Pa,Pb)- (2.2) 

a,b J ° J ° 

The PDF faui^at A*f) plays the role of a generalized number density to find a parton a 
carrying the momentum fraction x a of the parent hadron A with p a = x a p A denoting 
the four-momentum of the incoming parton to the hard scattering. We work in the 
QCD-improved parton model using the five-flavour scheme with Nf = 5 massless quarks 
q = u, d, c, s, b. The partonic subprocesses that contribute to the above scattering reaction 
at NLO can be generically written as 

a{Pa) + Kp b ) -> c{k c ) + d{k d ) {+e{k e )) , (2.3) 

where a, b, c,d,e G {g, u, d, c, s, b} if only weak corrections are considered. For photonic 
corrections also the photon has to be included as a possible external particle state. The 
additional emission of parton e appears in the real NLO correction to this process. The 
momenta k\ and k 2 in Eq. (12.11) emerge from the recombination procedure of the jet 
algorithm. At LO, however, we simply have k\ = k c , k 2 = ka, with fcx,i = kr s 2- Owing 
to the mass degeneracy (m q = 0) of the external quarks and the unitarity of the CKM 
matrix, the effect of the non-diagonal CKM structure vanishes in most contributions after 
taking the flavour sums. The only exception is the case of a W-boson exchange in the 
s-channel, where the different weights from the PDFs spoils the cancellation. However, 
this dependence turns out to be negligible and we set the CKM matrix to unity in our 
calculation. 
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The partonic subprocesses can be classified as follows: 



g + g -»■ 9 + 9, 




(2.4a) 


g + g^q + q, 




(2.4b) 


Ui + di — > Uj + dj, 




(2.4c) 


Ui + di^- Ui + d~i, 




(2.4d) 


qi + qi-* q'j + q'j, 




(2.4e) 


q + q ->■ q + q, 




(2.4f) 



where i, j denote the generation indices. We note that the processes listed in Eq. f)2.4p 
are merely representatives of a class that also include the reactions that are related by 
crossing symmetry. First, we categorize the processes according to the number of gluons 
and quarks that appear as external particles into the four-gluon (a), the two-gluon-two- 
quark (b), and the four-quark (c-f) processes. The four-quark processes can be subdivided 
into processes that involve a W-boson exchange diagram (c,d) and those that only contain 
neutral-boson exchange diagrams (e,f). A further distinction is made by distinguishing 
processes that involve both s-channel and t-channel diagrams (c,e), and those that only 
include either one (d,f). The LO Feynman graphs to the above process classes are shown 
in Fig. [TJ Furthermore, we can exploit the symmetry of the matrix element with respect to 
the interchange of the generation index of the first two quark generations. This reduces the 
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number of independent amplitudes that need to be evaluated and speeds up the numerical 
evaluation. 

The electroweak coupling constant is derived from the Fermi constant in the scheme 
via the following relation 



V2 



7T 



1 



(2.5) 



This input-parameter scheme avoids large logarithms of the light fermion masses generated 
by the running of the coupling constant a(Q) from the Thomson limit (Q = 0) to the 
electroweak scale (Q ~ Mw) and furthermore absorbs universal corrections induced by 
the p parameter (see e.g. Ref. (3l|). 

In order to describe the resonances of the intermediate vector bosons Z and W, we 
employ the complex-mass scheme [32|, 
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which fully respects gauge invariance. In this 
approach the square of the gauge-boson mass is defined as the position of the pole in 
the complex k 2 plane of the respective propagator with momentum k. The consistent 
replacement of the (squared) gauge-boson masses by complex values, 



M 2 v ^p\ = Ml- iM v T v , V = W, Z, 



(2.6) 



induces the adaption of all real quantities, 
mixing angle 8 W : 



In particular, this results in a complex weak- 



cos 



A^W • 2 a 

— 7T, sm 9 V 



(2.7) 



In order to ensure the correctness of the presented results two independent calcu- 
lations have been performed, resulting in two separate implementations for the numer- 
ical evaluation. Both calculations employ the Feynman-diagrammatic approach in the 
't Hooft-Feynman gauge for the loops and the Catani-Seymour dipole subtraction ap- 
proach 3J, |35[ to isolate and cancel infrared (IR) divergences. The results of the two 
calculations are in mutual agreement. 

In the first calculation all tree-level amplitudes are calculated and implemented by 
hand using the Weyl-van-der-Waerden spinor formalism as worked out in Ref. 36]. The 



virtual corrections are calculated using the Mathematica Packages FeynArts 3.6 [37 
and FormCalc 6.2 38]. The one-loop integrals are evaluated using a modified version 



of the LoopTools 2.4 [38| library, which was supplemented by the loop integrals with 
dimensionally regularized IR divergences that were not included in version 2.4. Addition- 
ally, an interface is implemented between LoopTools and the Collier library, which is 



based on the results of Refs. [39|, |40J for tensor and scalar one-loop integrals, respectively. 
This allows to utilize the unmodified code generated by FormCalc, while resorting to 
Collier for the evaluation of the loop integrals. The results of the two approaches are 
in perfect mutual agreement. The numerical integration is performed using the adapt- 
ive Monte Carlo algorithm Vegas 
chosen. 



411 . where a specific phase-space parametrization is 



The second calculation uses FeynArts 1.0 [42[] for generating the tree-level and one- 
loop diagrams and in-house Mathematica routines to obtain an analytic result which 
then is exported as Fortran source code. The loop integrals are evaluated using the 
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Collier loop library. The finite Catani-Seymour dipole subtraction terms and real 



emission matrix elements are built around amplitudes generated with the O'Mega [43 



matrix-element generator. Adaptive single-channel Monte Carlo integration is implemen- 



ted using the Vamp |44j library. Contrary to the first calculation, the second calculation 
does not implement subprocesses involving external bottom quarks which only amount 
to w 3% of the LO cross section and thus can be safely neglected at order a^a s . For all 
other subprocesses, both calculations are in excellent agreement for both the integrated 
and differential cross sections. 

2.2 Structure of the NLO calculation 

The calculation of the LO cross section is based on the full SM, i.e. all vector bosons, 
including the photon 7, are included. The Born diagrams for each process class defined 
in the previous section are shown in Fig. HJ 

At NLO we classify the corrections into photonic and purely weak corrections. This is 
possible, because each diagram in the full NLO correction, defined by the order a^a, con- 
tains exactly one electroweak gauge boson due to the single power in a. We can therefore 
uniquely assign each NLO contribution to either a photonic or a purely weak correction. 
The photonic corrections would constitute the O (a^a) corrections in a hypothetical gauge 
theory with the group S77(3)c x {7(1)qed- Therefore, they form a gauge-invariant subset 
and consequently so do the remaining purely weak corrections. In the following, we will 
use the notation O (a^a w ) to refer to the purely weak contributions. 

This gauge-invariant classification of the corrections allows us to tackle each class 
successively. In this work we present the calculation of the O (a^a w ) corrections, as 
motivated in the introduction. The inclusion of the photonic corrections in order to 
obtain the full EW corrections at the order oPa is left to the future. 



It is known [28|, |45| that partial cancellations can occur between the virtual weak 
corrections and the real emission of massive gauge bosons. However, to which extent this 
compensation takes place strongly depends on the experimental setup to reconstruct W 
and Z bosons. Considering that the weak-boson emission is a tree-level process which can 
be easily simulated with fully LO automatic tools, they are not further considered here. 

2.2.1 Virtual corrections 

The virtual corrections consist of the one-loop diagrams and the corresponding coun- 
terterms. Because we are restricting our NLO calculation to the order 0(a^a w ), only 
the interference terms shown in Fig. [2] are relevant. The generic diagram for the virtual 
corrections to the process class (12.4bj) is shown in Fig. [2](a). The corrections constitute 



the purely weak O (a w ) correction to the LO O (0%) cross section; some representative 
diagrams are depicted in Fig. [2£a'). In case of the process classes f)2.4d -f) with four ex- 
ternal quarks, however, we have LO amplitudes of the order a w and a 8 . This leads to the 
two types of virtual corrections shown in Figs. [2](b,c). Here, we can identify genuine QCD 
corrections, such as the contributions shown in Fig. Et^c') and the first vertex correction 
in Fig. EJ^b'), and weak corrections, such as the second vertex correction in Fig. EJ^b'). 
However, there are also contributions such as the box diagram in Fig. [2](b') that cannot 
be assigned uniquely to QCD nor to weak corrections. This indicates that the separation 
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Figure 2: The virtual corrections of O (a^a w ) illustrated in terms of interference diagrams 
of generic Feynman graphs (a-c) and a set of sample diagrams below (a'-c'). The white 
circles and the double-circles in the interference diagrams represent tree-level and one-loop 
subgraphs, respectively. 
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Figure 3: One-loop and counterterm diagrams for the process class (12.4b|) grouped into 
the self-energy (a), vertex (b), and box (c) corrections. All counterterms are restricted to 
the order O (a w ). V denotes the vector bosons W and Z, and q' the weak-isospin partner 
of q for V = W and q' = q for V = Z. 



of the "QCD corrections" from the "weak corrections" is not properly defined and instead, 
one must treat them together as a whole, defined by the order in perturbation theory. The 
contributions we referred to by the "QCD corrections" above contain infrared divergences 
and, therefore, the real emission of an additional gluon must be considered, which will be 
discussed in the next section. A more complete set of one-loop and counterterm diagrams 
associated with the schematic illustrations in Figs. [2(a), [2(b), and |2](c) can be found in 
Figs. [3j [U and [5j respect ivelyfj] 

Ultraviolet (UV) divergences are regularized dimensionally. For the IR singularities 
also dimensional regularization is used, but our second calculation optionally employs 



infinitesimal masses as regulators. It has been shown in the appendix of Ref. [46[ that 
rational terms of IR origin cancel in any unrenormalized scattering amplitude, so that they 
need not be further considered in the calculation of the one-loop amplitudes. The only 
remaining source of rational terms of IR origin involve the wave-function renormalization 
constants, which were calculated separately and are given below. 

The external fields are renormalized in the on-shell scheme and hence, all self-energy 
corrections to the external (on-shell) legs vanish and can be omitted already at the level 
of diagram generation. The renormalization of the strong coupling constant is done using 
the MS scheme for the Nf = 5 light quarks and by subtracting the contribution of the 



1 In the case with external bottom quarks there exist additional diagrams in Fig. [3] and Fig. QJb) 
due to the no n- vanishing mass of the weak-isospin partner t, where in place of a charged would-be 
Goldstone boson ^ is exchanged. 
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Figure 4: One-loop and counterterm diagrams for the processes ( 12.41) (d-f) of O (a s a w ) 
grouped into vertex (a,b) and box (c) corrections. The triangle insertions are further 
subdivided into QCD (a) and weak (b) corrections, and consequently the associated coun- 
terterms are restricted to the order O (a s ) and O (a w ), respectively. V denotes the vector 
bosons W and Z. 
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Figure 5: One-loop and counterterm diagrams for the processes (I2.4p (d— f) of O (cCg) 
grouped into self-energy (a), vertex (b), and box (c) corrections. All counterterms are 
restricted to the order O (a s ). 



S 



heavy top-quark loop in the gluon self-energy at zero momentum transfer. Therefore, the 
running of the strong coupling constant is driven by the five light quark flavours only. 

Similar to the loop diagrams, the renormalization constants need to be evaluated at 
different orders, where the ones of O (a s ) are 



SZ ga \o(a s ) 
SZ q \ {a s ) 



Air 



T"t) (A- + ln(g))-iB„(0, mt , mt ) 



2ir 



|-^r) 5 (o,o,o)-ifi c 



3tt 



B (0, 0,0), 



and the ones of O (a w ) are 



5Zg\ (a„) 



a 

Air 

a 

Air 



(9a 



2£? 1 (0,0,/i z )) 



2B 1 (0,0,/z w )) 



(g+) 2 (1 + 2^(0,0,^)) 



(2.8) 

(2.9) 
(2.10) 

(2.11) 
(2.12) 



where our notation for the 2-point functions B 01 follows Ref. [39|, // is the arbitrary 
reference mass of dimensional regularization, the renormalization scale, and 



A, 



A-D 



7e + ln(47r) 



(2.13) 



r> it 

denotes the standard one- loop UV divergence in D dimensions. Here 5Z q and SZq 
are the field- renormalization constants of the right /left-handed quark fields and of the 
gluon field, respectively, and 5Z gs connects the bare (g s ,o) and the renormalized (g s ) strong 
coupling constant, g s>0 = (l+SZ gs )g 8 . The couplings g^ are defined via the 3rd component 
of weak-isospin, I^ q , and the electric charge Q q of the quark q, 



I 3 

w,q 



slQ q 



-1/2, +2/3), 
-1/2,-1/3), 



q 



di. 



(2.14) 



At this point some comments on the use of complex masses and couplings, as dictated 



by the complex- mass scheme 32|, [33[, are appropriate. This scheme, which was primarily 
introduced to achieve a consistent, gauge- invariant description of gauge-boson resonances 
at LO and NLO, does not only comprise the consistent use of complex parameters in amp- 
litudes, but also complex generalizations of the renormalization constants for the complex 
masses and couplings as compared to on-shell renormalization schemes for real masses (see 
e.g. Ref. 47J). Note, however, that the order O (c^a) of the corrections considered in this 
calculation does not involve weak corrections to gauge-boson propagators and weak coup- 
lings, so that the complex generalization of the relevant renormalization constants given 
in Eqs. ( I2.8j) — (12. 12[) just concerns the insertion of complex masses and weak couplings. 



2.2.2 Real corrections 

The real corrections receive contributions from the partonic processes that are illus- 
trated in terms of interference diagrams in Fig. El The contribution shown in Fig. E{a) 
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(a) gluon emission 



(b) gluon induced 



Figure 6: The real corrections at O (a^a w ) illustrated in terms of interference diagrams 
of generic Feynman graphs. The circles represent tree-level subgraphs. 



corresponds to the additional emission of a gluon from the processes (12.4H — f). properly 
taking into account only the interference terms that contribute at the order 0(a^a w ). 
The gluon-induced corrections in Fig. Mb) are obtained from the preceding by crossing 
the gluon into the initial state. 

The real-emission cross section contains infrared divergences in the phase-space in- 
tegration which have their origin in the regions where a final-state parton becomes soft 
or collinear to another parton. The soft and the final-state collinear singularities cancel 
against the corresponding singularities in the virtual corrections for sufficiently inclusive 



observables by virtue of the Kinoshita-Lee-Nauenberg theorem [48|, |49(. The remaining 
initial-state collinear singularities are process independent and absorbed into the NLO 
PDFs by QCD factorization, which is technically accomplished by subtracting a so-called 
collinear counterterm (da ) from the NLO cross section. The subtraction formalism re- 
shuffles the IR singularities by constructing a subtraction term (da A ) to the real correction 
(dcr R ) which mimics its singular behaviour to render their difference integrable in four di- 
mensions. The subtraction term can be integrated analytically in D = 4 — 2e dimensions 
over the singular one-particle subspace, generating 1/e and 1/e 2 poles that cancel against 
the corresponding poles in the virtual corrections (da v ) and the collinear counterterm. 



To this end, we employ the Catani-Seymour dipole subtraction formalism [34|, which 
constructs the subtraction term in terms of so-called dipoles, which are built from the 
LO amplitudes (<r B ) and dipole operators (dVdi po ie) which in general contain colour and 
helicity correlations. The NLO contribution to the hard-scattering cross section from 
Eq. (12.21) can be schematically written as, 



a 



NLO 



da 1 



da N + da 

J2 



C 



da v + da° + / da 



/[(^ R ) £=0 -(d- A ) e= 

J ( d(jR ) e =o - ( dcjB ® dv ^ 

f dx [ da B g) (K(x) + P(x,f4)) 

J2 



dipolcs 



dcr v + da a (8) /(e) 



(2.15) 
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0(a s ) 



Figure 7: The building blocks for the subtraction terms at O (a^a w ) illustrated in 
terms of interference diagrams of generic Feynman graphs. The circles represent tree- 
level subgraphs. The symbol (g) denotes possible colour and helicity correlations between 
the dipole operators and the LO amplitudes. 



where f denotes the integration over the m-particle phase space, (g) encodes the possible 
colour and helicity correlations between dcr B and the dipole operators, and the integration 
over x corresponds to a convolution over the momentum fraction of the incoming partons. 
The insertion operators /(e), K(x), and P(x,fip) emerge from the collinear counterterm 
and the integration of the subtraction term over the singular one-particle subspace. More 



details and the explicit expressions for /(e), K(x), and P(x, /i|) can be found in Ref. 34 
The dipole operators introduced above are all of 0(a s ), which in turn forces us to 

restrict the calculation of the colour-correlated LO amplitudes to the interference terms 

of O (a s a w ). This is schematically illustrated in Fig. [7J 

For the second calculation which used both dimensional and mass regularization for 

the IR divergences, the generalization of the above formalism for massive partons 35[ is 

used for the latter scheme. 

3 Numerical results 

3.1 Input parameters and setup 



For the numerical evaluation we use the input parameters of Ref. [50 



= 1.16637 x 1(T 5 GeV" 2 , a s (M z ) = 0.129783, 

M° s = 80.398 GeV, M° s = 91.1876 GeV, 

= 2.141 GeV, r£ s = 2.4952 GeV, 

m t = 172.5 GeV. (3.1) 

Note that there is no Higgs-boson mass dependence in the predictions for dijet production 
at the considered order. The complex-mass scheme for the W and Z bosons corresponds 
to a fixed-width description of the resonance and requires a conversion of the on-shell 



gauge-boson masses to the pole masses My as follows 51 



cv ' " cy ' v V \ M v S J 




M os r os 
Mv = ^, rv = ^, c y = Wl+(-^) , (3.2) 
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giving 



M w = 80.3695 . . . GeV, M z = 91.1535 . . . GeV, 

r w = 2.1402... GeV, T z = 2.4943 ... GeV. (3.3) 



For the PDFs we use the CTEQ6L1 [52J set, which dictates the value of a s (Mz) in 
Eq. (13. ip . A consistent QCD calculation, of course, requires the use of LO, NLO, or 
NNLO PDFs for the calculation of the respective cross-section prediction. Our aim here 
is, however, to provide a relative correction factor for weak effects that is to be applied to 
state-of-the-art QCD predictions. This factorization procedure is better motivated than 
just adding QCD and electroweak corrections because of the factorization of IR-sensitive 
corrections. The weak correction factor has to be derived from a single PDF set, but is 
rather insensitive to PDFsJ§ The renormalization scale is set equal to the factorization 
scale, chosen as the transverse momentum of the leading jet 

A*R = A*F = \i = fc T ,i. ( 3 - 4 ) 
3.2 Phase-space cuts and event selection 

The definition of an IR-safe jet observable requires the recombination of soft and/or 
collinear partons in the final state, as well as constraining the phase space by imposing 



cuts. The jets emerge from the final-state partons via the anti-A>p algorithm |53j, where 
we have set the angular separation parameter to R — 0.6. Recombination is performed 
using four-momentum summation. 

We require the jets to have transverse momenta larger than k^ et and demand 
them to be central by restricting their rapidities yi to the range \yi\ < y^t with the values 

y$ = 2.5, kSft* = 25 GeV. (3.5) 

3.3 Results 

In the following, we present the numerical results for dijet production at the LHC, i.e. 
for a pp initial state, at the centre-of-mass (CM) energies of y/s = 7 TeV, 8 TeV, and 
14 TeV, and at the Tevatron, i.e. for pp collisions at the CM energy of y/s = 1.96 TeV. 

We denote the full LO cross section through O ( 2 ) by a and define the NLO 

corrections relative to the LO cross section via 

^^"xll + O- (3-6) 

Furthermore, we denote the LO QCD cross section of O (a 2 ) by 0q CD and introduce a 
correction factor for the left-over LO contributions of order a s a and a 2 , 

a° = a° QCD x(l + 5^). (3.7) 



2 The consistent use of LO PDFs is fully justified for the electroweak tree-level contributions and 
the weak loop corrections (containing the large weak Sudakov logarithms) to the QCD channels. The 
remaining part of the calculated correction actually has the character of a QCD correction, viz. the inter- 
ference of QCD loops with weak LO diagrams and the corresponding real corrections, strictly requiring 
NLO PDFs. Recall, however, that the two different loop contributions cannot be separated, as explained 
in Sect. 12.21 Since the calculated corrections are generally rather moderate, our procedure is certainly 
acceptable within the remaining uncertainty of the aimed combination of QCD and weak corrections. 
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PP — > 3 3 + X at y/s = 7 TeV 



A/12/Gev 


50-oo lOO-oo 200 - 00 500 - 00 lOOO-oo 2000 - 00 3000 - 00 


cr /nb 

CT QCD/ nb 


78600(2) 25496(1) 3879.2(3) 80.807(2) 2.49314(5) 33.0487(9) • 10~ 3 890.83(1) ■ 10" 6 
78561(2) 25490(1) 3878.07(10) 80.650(3) 2.4756(2) 32.4826(8) • 10~ 3 870.69(2) • 10" 6 


<5ew/% 

E<V% 


-0.02 -0.03 -0.09 -0.39 -1.09 -2.48 -3.60 
0.05 0.01 0.04 0.20 0.70 1.75 2.31 
0.02 -0.02 -0.05 -0.19 -0.38 -0.73 -1.29 


<7 NLO /nb 


78581(2) 25487(1) 3875.7(3) 80.493(2) 2.46627(6) 32.244(1) ■ 10^ 3 859.48(4) ■ 10" 6 
PP — ► 33 + X at s/s = 8 TeV 


Mi 2 /GeV 


50 - 00 IOO-00 200 - 00 500 - 00 1000 - 00 2000 - 00 3000 - 00 


cr /nb 


94683(2) 31201(1) 4868.2(3) 108.331(7) 3.65684(7) 61.498(1) • 10" 3 2.45475(6) • 10" 3 
94638(2) 31197(2) 4867.1(1) 108.152(3) 3.6344(1) 60.5237(9) • 10~ 3 2.40056(5) • 10~ 3 


el-loop /<w 
d weak / /o 

E<V% 


-0.02 -0.03 -0.08 -0.37 -1.05 -2.44 -3.62 
0.04 0.01 0.03 0.17 0.61 1.61 2.25 
0.02 -0.02 -0.05 -0.20 -0.44 -0.83 -1.37 


(7 NLO /nb 


94661(2) 31191(1) 4864.1(3) 107.926(8) 3.61864(10) 60.019(2) • 10" 3 2.3678(1) • 10" 3 
PP — ► 33 + X at y/s = 14 TeV 


Mi 2 /GcV 


50 - 00 IOO-00 200 - 00 500 - 00 IOOO-00 2000 - 00 5000 - 00 


cr /nb 
°QCD/ nb 


198483(4) 69335(3) 11858.2(8) 334.17(1) 14.9435(4) 456.64(1) • 10~ 3 954.71(2) • 10" 6 
198410(5) 69329(2) 11856.0(10) 333.849(8) 14.8938(3) 452.12(2) • 10 -3 930.90(3) ■ 10~ 6 


O op /% 
E<V% 


-0.02 -0.03 -0.07 -0.31 -0.88 -2.20 -5.53 
0.03 0.01 0.02 0.10 0.34 1.00 2.56 
0.01 -0.02 -0.05 -0.22 -0.55 -1.20 -2.98 


(7 NLO /nb 


198444(4) 69315(3) 11849.5(8) 333.12(1) 14.8117(4) 446.69(2) • 10~ 3 903.19(7) • 10" 6 



Table 1: Integrated dijet cross sections and respective corrections for various ranges of 
the dijet invariant mass Myi at the LHC with CM energies 7 TeV, 8 TeV, and 14 TeV. 

For the NLO cross section this leads to 

^=<^x(i+^)x(i+*i5 p ) 

^4dx(1+« + C P )- (3-8) 
With respect to the LO QCD cross section, the total correction is given by the sum 
^wcak P + ^ew- Owing to the rather moderate size of the corrections, the difference in 
defining the NLO corrections ^ak 13 relative to a or ctq CD constitutes a higher-order 
effect which is negligible. 

3.3.1 The dijet invariant mass at the LHC 

The dijet invariant mass is defined as M12 = \J {E\ + E 2 ) 2 — (pi + P2) 2 , where E\ )2 
and pi t 2 denote the energies and the momenta of the leading and subleading jets, respect- 
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ively. In Table [T] we present the integrated LO cross sections cr", c^qcd; the NLO cross 

section cr NLO , and the relative correction factors <5^ e 1 °£ p , 5ew> Yl $ = ^EW+^woak P f° r differ- 
ent cuts on the dijet invariant mass M 12 at the LHC. We note that the lowest cut value 
of Myi > 50 GeV is already covered by the standard setup via the cut on the transverse 
momenta of the jets, k^ ct = 25 GeV, and does not constitute a further restriction on the 
cross section. 

The cross section is dominated by the region close to the cut, which is reflected by the 
rapid decrease of the integrated cross section with increasing values for the cut on Mi 2 . 
Therefore, also the corrections are dominated by the region given by the lowest accepted 
M 12 values. 

Comparing the two cross sections <t° and 0q CD , we find that the LO cross section with 
minimal cuts is predominantly given by the QCD cross section and that the electroweak 
effects (5ew) typically stay below the per-cent level. However, we observe a steady increase 
in with higher M 12 , which can be explained by the parton luminosities at the LHC: 
As can be seen from Fig. [SJ the LO cross section is dominated at lower M 12 by the gg- and 
gq- initiated processes, for which 5^ vanishes. At higher M 12 , the gg-initiated processes 
with (JJlyy ^ become dominant, leading to the behaviour described above. The running 
of the strong coupling also acts in favour of increasing <5 E r {y with higher cuts on Mi 2 . 
Comparing the distributions in Fig. [8] for y/s = 7 TeV and 14 TeV we observe a shift of 
the transition region from gg- to gg-domination from ~ 1.5 TeV to ~ 3 TeV, respectively. 
This trend is due to the fact that lower values of y/s require larger partonic momentum 
fractions x for a fixed value of Mu, and the (valence) quark PDFs are enhanced over the 
gluon PDF for larger x. This also explains the y/s dependence of the EW contribution 
^EW' w hich decreases with higher CM energies for the same cut on M 12 . 

The purely weak corrections are negative throughout and increase in magnitude from 
-0.02% to -3.6% in case of the = 7 TeV setup for a M 12 cut of 50 GeV and 3 TeV, 
respectively. This behaviour partly originates from the corrections containing weak logar- 

/ M 2 \ 

ithms In ( j , which become larger by effectively restricting the cross section to higher 
energy scales Q via increasing cuts on the invariant mass. 

Compared to the relative weak corrections &^!°£ p show a far weaker dependence 
on the CM energy. The corrections 5^ and (^!°£ p are similar in magnitude, but of 
opposite sign, leading to large cancellations in the sum. 

Note that the weak loop corrections <^!°£ p in the TeV range are small in comparison 

2 ( M 2 \ 

to the typical size of the EW Sudakov factor ^ In ( j , which amounts to tens of per 

cent. This is due to the fact that also for large invariant mass M 12 the cross section is 
not dominated by the Sudakov regime which requires that the absolute values of both 
partonic Mandelstam variables s = (p a + Pb) 2 and t = (p a — k c ) 2 are much larger than 
M w . Instead, the cross section here is dominated by the Regge (forward) region where s 
is large, but \t\ remains small. 

This feature is also evident in the dijet invariant-mass distributions shown in Fig. |9j 
The relative corrections given in Table [1] are almost identical to the corrections in Fig. [9] 
at the respective cut value of M 12 , since the corrections are dominated by the region close 
to the cut. 
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M12 [GeV] M12 [GeV] 

Figure 8: The LO contributions to the dijet invariant mass M 12 distribution from the 
different initial-state parton combinations at the LHC with CM energies 7 TeV, 8 TeV, 
and 14 TeV. Left: absolute predictions; right: relative contributions S tree . 
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da°/dM 12 [nb/GeV] 6 [%] 
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Figure 9: Differential distributions with respect to the dijet invariant mass M\i at 
LHC with CM energies 7 TeV, 8 TeV, and 14 TeV. Left: absolute predictions; ri; 
relative contributions 5. 
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d 2 cr°/(dM 12 dy*) [nb/GeV] 
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Figure 10: Double-differential distribution with respect to the dijet invariant mass M 12 
and y* at the LHC with a CM energy of 8 TeV. In the absolute prediction (uppermost 
plot) the cross section is divided by the bin width in y*. 
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Figure 11: Same as in Fig. [101 but for a CM energy of 14 TeV. 
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Figures [TOl and [TTI show the results for the Myi distributions corresponding to different 
rapidity ranges. Specifically, the phase space is divided into different regions of y*, which 
is defined as half the rapidity difference of the two leading jets, 

✓ = (3.9) 

For the binning in y* we have chosen five bins of equal width covering the range < y* < 
2.5. In all y* bins we observe the above described behaviour of the relative corrections that 
increase in magnitude with higher invariant mass, where the tree-level EW contributions 
and the radiative weak corrections give positive and negative contributions, respectively. 
Furthermore, a strong dependence on y* is observed, such that both, 5^ and 5^?°^, are 
larger for smaller values of y*. For y/s = 8 TeV and the invariant mass of 3 TeV the 
corrections 5^ and ^!°£ p in the first bin (y* < 0.5) amount to approximately 15% and 
—10%, respectively, whereas in the highest y*-bin (2 < y* < 2.5) we observe about 1.5% 
and —2.5% for the respective corrections. As it is evident from the numbers quoted above, 
decreases more rapidly than <5^ C 1 °£ P for higher y*. The sum of both contributions, 
thus, has a positive net-contribution for y* < 0.5, gradually decreasing with higher y* 
values. In the range of highest y* (2 < y* < 2.5), even a negative net-contribution results. 
Qualitatively we observe the same behaviour also at the CM energy of \/s = 14 TeV. Here 
we examine dijet invariant masses up to 6 TeV, where we observe approximately 20% and 
-15% in the first bin (y* < 0.5), and 2% and -5% in the highest y*-bin (2 < y* < 2.5) 
for 6%fi and ^°° p , respectively. 

As already indicated in the discussion of integrated cross sections above, the behaviour 
of the corrections can be better understood by identifying the Sudakov regime in phase 
space. For the contributions with 2 — s- 2 kinematics, i.e. the LO cross sections and the 
virtual corrections, the following relations hold: 

\m\ = \-m\ = y\ § = m 2 12 , t = ~ TT ^ F , ^-^--L (3.io) 

where t/i and s, t, u denote the rapidity and the partonic Mandelstam variables. The 
two different signs refer to y\ ^ 0. For large M 12 3> Myv an d small y* we reside in the 
Sudakov regime, where all scales are simultaneously large compared to the vector-boson 
masses. Here, all logarithms of the form a w ln 2 (Q 2 /My), with Q 2 = s, — i, —u become 
large, leading to the observed enhancements of the weak radiative corrections. As can be 
seen in the absolute distributions in Figs. [TU] and [TTJ however, this regime delivers only 
a small fraction to the cross section if the high-energy region is defined via cuts on the 
invariant mass M\2- 

Although the corrections snow a similar behaviour, their origin is of a completely 
different nature: As discussed above, the bulk of the contributions to 5^ originate from 
the gg-initiated processes, in particular from the partonic subprocesses uu — > uu and 
ud — > ud. Owing to their colour structure, the interference terms of O (a s a) considered 
here receive contributions from products of t-, u-, and s-channel diagrams only, but not 
from squares of those topologies. Therefore, these interferences appear to be more central 
than the larger LO QCD contributions which are dominated by squared topologies. 

In summary, for high M 12 , where the corrections are largest, the cross section is dom- 
inated by the contribution coming from the highest possible y* values. However, this 
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PP — > 3 3 + X at y/s = 7 TeV 



fc T ,i/GcV 


25 -oo 50-oo 100 -oo 200 - oo 500 - oo 1000 - oo 1500 - oo 


cr /nb 
^qcd/ 1 * 


78600(2) 5417.1(1) 291.205(9) 11.1765(2) 63.697(2) ■ 10~ 3 396.843(8) • lO" 6 5.9746(10) ■ 10~ 6 
78561(2) 5413.8(2) 290.517(9) 11.0873(3) 61.251(2) • 10~ 3 353.47(1) • 10" 6 4.9611(9) • 10~ 6 


rl-loop 10/ 
d weak / /o 

EV% 


-0.02 -0.10 -0.34 -0.99 -2.96 -5.12 -6.13 
0.05 0.06 0.24 0.80 3.99 12.28 20.44 
0.02 -0.04 -0.10 -0.19 1.04 7.16 14.31 


a NLO /nb 


78581(2) 5411.9(1) 290.216(9) 11.0664(3) 61.887(4) ■ 10" 3 378.76(4) ■ 10~ 6 5.6705(8) • 10~ 6 

PP — ► 33 + X at y/s = 8 TeV 


fc T ,i/GeV 


25 -oo 50-oo 100 -oo 200 - oo 500 - oo 1000 - oo 1500 - oo 


cr°/ n b 
°QCD/ nb 


94683(2) 6762.3(2) 380.68(1) 15.7387(5) 105.965(3) ■ 10~ 3 877.88(2) ■ 10" 6 20.4538(3) ■ 10~ 6 
94638(2) 6759.0(3) 379.90(1) 15.6330(6) 102.626(2) • 10~ 3 795.70(3) • 10~ 6 17.3482(2) • 10~ 6 


Sew/% 

EV% 


-0.02 -0.09 -0.33 -0.97 -2.96 -5.35 -6.70 
0.04 0.05 0.21 0.68 3.26 10.33 17.90 
0.02 -0.04 -0.12 -0.29 0.30 4.98 11.20 


a NLO /nb 


94661(2) 6756.1(2) 379.44(1) 15.5871(6) 102.932(5) ■ 10~ 3 835.29(9) ■ 10~ 6 19.292(3) • 10" 6 
PP — ► 33 + X at y/s = 14 TeV 


fc T ,i/GeV 


25 -oo 50-oo 100 -oo 200 - 00 500 - oo 1000 - oo 2500 - oo 


cr°/ n b 
^QCD/nb 


198483(4) 16194.3(5) 1074.11(3) 56.405(2) 671.10(2) ■ 10~ 3 12.0383(4) ■ 10~ 3 8.8504(2) • 10~ 6 
198410(5) 16189.7(5) 1072.85(4) 56.204(1) 661.52(2) ■ 10 J 11.5060(2) ■ 10 J 7.4826(2) • 10 


E<V% 


-0.02 -0.08 -0.28 -0.84 -2.72 -5.48 -10.49 
0.03 0.03 0.12 0.36 1.44 4.62 18.28 
0.01 -0.05 -0.15 -0.48 -1.28 -0.86 7.79 


a NLO /nb 


198443(4) 16181.6(5) 1071.15(3) 55.935(2) 653.10(3) • 10~ 3 11.4076(10) • 10~ 3 8.066(1) • 10~ 6 



Table 2: Integrated dijet cross sections and respective corrections for various ranges of 
the transverse momentum of the leading jet kr,i at the LHC with CM energies 7 TeV, 
8 TeV, and 14 TeV. 

region receives the smallest corrections, leading to the small corrections that we observe 
in Table [1] even for very high M\2 cuts. As we will see in the next section, this behaviour 
is reversed in case of the transverse-momentum distributions. Despite the different origin 
of the relative corrections 5^l°^ p and 5gw: they conspire together to large cancellations in 
the sum ^EW+^weak P - The degree of this compensation, however, depends on the chosen 
cuts defining the observable. 

3.3.2 The transverse momenta of the leading and subleading jet at the LHC 

In Table [2] we list the various integrated cross sections at LO, the NLO cross section, 
and the correction factors defined above for different cuts on the transverse momentum of 
the leading jet, k T ^. The cut fcx.i > 25 GeV in the first column is already imposed by the 
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default set of cuts and does not represent a further restriction to the cross section. Again, 
the integrated cross section decreases rapidly with more restrictive cuts on the transverse 
momentum, indicating that the cross section, and with it the corrections, are dominated 
by the region with the lowest accepted /ct,i- 

The weak radiative corrections display the expected Sudakov-type behaviour with 
increasing negative corrections for higher fcx.i-cuts and only a modest dependence on the 
CM energy of the collider. For cut values of 25 GeV up to 1.5 TeV they increase from 
-0.02% to -6% for yfs = 7 TeV. For y/s = 14 TeV and a cut of fc T) i > 2.5 TeV the 
radiative corrections even amount to —11%. The LO EW contributions show a much 
stronger dependence on the collider energy than the loop corrections, in particular for 
more restrictive cuts. This \/s dependence is also stronger than the one observed for the 
dijet invariant- mass spectra discussed in the previous section. 

Considering that &t,i < M 12 /2 at LO and that the cross section as well as the correc- 
tions are dominated by the region of lowest accepted M 12 , one might naively expect that 
the results for a fixed cut on fcx,i should be comparable to the respective corrections for 
a M 12 cut of twice that value. However, the region close to M 12 ~ 2A>r,i at LO requires 
central jet production in the partonic CM frame, i.e. y* ~ 0, so that the cross section 
defined via the /ct,i cut is dominated by the Sudakov regime. This explains the smaller 
cross section and the larger corrections as compared to the cross section defined via the 
corresponding M 12 cut discussed in the previous section. 

The differential distributions in &>r,i shown in Fig. [12] cover the range up to &t,i = 
1.5 TeV for \/s = 7 TeV, 8 TeV, and up to /ct,i = 3 TeV for y/s = 14 TeV and further 
underline the observations made above. The weak corrections display the expected beha- 
viour from the Sudakov logarithms with corrections reaching up to —6% at /ct,i = 1-5 TeV 
for yfs = 7 TeV, 8 TeV, and -12% at fc T ,i = 3 TeV for y/s = U TeV. On the other hand, 
the tree-level corrections increase with higher fcx,i and reach approximately 16% at 
fc T ,i = 1.5 TeV for ^ = 7 TeV, 8 TeV, and 20% at fc T ,i = 3 TeV for ^ = 14 TeV, 
resulting in significant cancellations in the sum ^EW+^weak P - 

By introducing a further binning in y* we obtain the double differential distributions 
shown in Figs. [13] and [Ml At higher values of the transverse momentum, the production 
of the jets is required to be more and more central in the partonic CM frame, leading 
to the observed rapid decrease in the cross section for the bins with higher values of y* . 
In contrast to the M 12 distribution, the bin with the smallest value for y* is the most 
dominant in the high-/cx,i tail. Moreover, both the tree-level EW corrections 5^ and 
the one-loop weak radiative corrections <5^e°£ p are only slightly affected by the y* binning, 
while there is a significant dependence in the M 12 distribution discussed in the previous 
section. 

Additionally, we present the corresponding transverse-momentum distributions with 
respect to the subleading jet in Figs. [[5] and [£7] Recall that leading and subleading jets 
have the same transverse momenta (/ct,i = ^t,2) i n all 2 — > 2 particle configurations, 
i.e. that only the real emission corrections to the four-quark channels show a different 
behaviour here. In particular, <5j^ remains the same when going from the leading to the 
subleading jet. On the other hand, the weak loop corrections <^g°£ p turn out to be more 
pronounced for small y*. This is due to the fact that subleading jets fill bins with smaller 
A/p in the spectra which rapidly decrease with higher k^s. Since the positive real-emission 
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Figure 12: Differential distributions with respect to the transverse momentum of the 
leading jet fcx,i at the LHC with CM energies 7 TeV, 8 TeV, and 14 TeV. Left: absolute 
predictions; right: relative contributions S. 
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Figure 13: Double-differential distribution with respect to the transverse momentum of 
the leading jet /ct,i and y* at the LHC with a CM energy of 8 TeV. In the absolute 
prediction (uppermost plot) the cross section is divided by the bin width in y*. 
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Figure 14: Same as in Fig. [131 but for a CM energy of 14 TeV. 
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Figure 15: Differential distributions with respect to the transverse momentum of the 
subleading jet /ct,2 at the LHC with CM energies 7 TeV, 8 TeV, and 14 TeV. Left: 
absolute predictions; right: relative contributions 5. 
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Figure 16: Double-differential distribution with respect to the transverse momentum of 
the subleading jet k^^ and y* at the LHC with a CM energy of 8 TeV. In the absolute 
prediction (uppermost plot) the cross section is divided by the bin width in y*. 
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Figure 17: Same as in Fig. [161 but for a CM energy of 14 TeV. 
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Figure 18: The LO contributions to the dijet invariant mass Mi 2 distribution from the 
different initial-state parton combinations at the Tevatron. Left: absolute predictions; 
right: relative contributions 5 trcc . 



contribution on the k^ axis of the subleading jets is, thus, shifted to the left as compared 
to the leading jet, the sum of negative virtual and positive real corrections for a fixed bin 
is somewhat shifted to more negative values for the subleading jet. This effect is largest 
for the smallest y* where the cross section is largest, because the real emission has a 
particular tendency to reduce maxima in distributions. 



3.3.3 Dijet production at the Tevatron 

In the following, we present the results for dijet production at the Tevatron, i.e. for a 
pp initial state and a CM energy of y/s = 1.96 TeV. 

Figure [18] shows the various contributions of the different partonic channels at LO 
contributing to dijet production at the Tevatron for the invariant-mass distribution of 
the two jets. Since the Tevatron is a pp collider, with valence quark-antiquark pairs in 
the initial state, there is a strong qq dominance at large values of the invariant mass M 12 
(several hundred GeV) , which requires large scattering energies and thus large momentum 
fraction x of the partons. At the moderate values M12 ^ 500 GeV, there is still some 
dominance of channels with gluons in the initial state, with even the largest contribution 
from gg scattering for very low M 12 , because the gluon PDF has the strongest rise at 
small x. 

The differential distributions with respect to the dijet invariant mass M12, the trans- 
verse momentum of the leading jet /ct,i and the subleading jet £>x,2 are shown in Fig. [T9Ta- 
c), respectively. In accordance with the observations made for the LHC, the weak correc- 
tions 5^!°k P are much smaller for Mi 2 -based observables as compared to those based on the 
transverse momenta of the jets. We further observe that the LO EW contributions 5^ in 
the M12 distribution is similar in magnitude and opposite in sign as compared to S^^ p , 
similarly to the case of the LHC. For the transverse-momentum distributions, however, 
becomes negative for higher /ct,i, further increasing the corrections in ^Ew+^wcak 13 ' 
which reach around —12% for k^^ = 800 GeV. It is interesting to note that the radiative 
corrections #weak P are similar for both, the leading and subleading jet /ct, which is different 
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Figure 19: Differential distributions with respect to the dijet invariant mass M\i (a), 
the transverse momentum of the leading jet /ct,i (b) and the subleading jet /ct,2 (c) at 
the Tevatron for a CM energy of 1.96 TeV. Left: absolute predictions; right: relative 
contributions S. 
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from the behaviour we observed at the LHC, where the corrections to &t,2 are signific- 
antly larger. To understand this, recall that the difference in the corrections for leading 
and subleading jets is merely due to real emission corrections in the four-quark channels. 
These are the real QCD corrections to the interference of weak and QCD tree diagrams. 
Generically these interferences, which are part of 5^ in LO, are much smaller at the 
Tevatron than at the LHC, as is obvious from 5jfw in Fig. [19j Moreover, it is interesting 
that the corrections are even slightly larger for the leading jet than for the subleading 
one, because the real QCD corrections to the weak-QCD interference are negative at high 
energy scales, as also observed for 5^ before. 

The M12 and /ct,i distributions with a further binning in y* are shown in Figs. [20] and 

EU 



3.3.4 Comparison to other work 



Preliminary results for the weak radiative corrections to dijet production at the LHC 
have been presented by A. Scharf in the proceedings contribution 18j], where the contri- 
butions from external bottom quarks were not considered as part of dijet production, but 



discussed separately. For comparison, we here adopt the calculational setup of Ref. [18 
for the LHC, 



= 14 TeV, 
y cut : none, 
a s = 0.1, 
M w = 80.425 GeV, 
PDF set : CTEQ6L, 



U cut 
/l T1 

A*f = Mr = 



- /i T 2 
ol,cut 

a 
M z 



50 GeV, 
100 GeV, 
1/128, 

91.1876 GeV, 



(3.11) 



which is partially inferred from Ref. [19|. According to the author of Ref. [18] . further 
details on the jet algorithm and the precise treatment of the W/Z resonances are not 
available anymore, but those loose ends should only play a minor role. This is confirmed 
in Fig.[22]which shows our result on the corrections to the transverse-momentum spectrum 
of the leading jet in comparison to the one shown to Fig. 9 of Ref. 18j. The two results 
show good agreement over the considered &t,i range. 



We have also tried to perform a tuned comparison to the results of Ref. [17| , but have 
not found agreement El 



4 Conclusions 



In this paper we have presented the calculation of the most important electroweak 
corrections to dijet production at the LHC and the Tevatron. These corrections com- 
prise electroweak contributions of O (a s a, a 2 ) to the LO QCD prediction as well as NLO 
corrections through the order a 2 a. Guided by the electroweak Sudakov-type logarithms 
induced by soft or collinear W/Z exchange at high energies, we have restricted ourselves 
to the calculation of the purely weak loop corrections in a first step. 

3 To find the source of discrepancy seems to require a careful comparison of individual components of 
the calculation. Since the correctness of our results is backed by our two calculations and the comparison 
to Ref. [3] , we do not see a reason to await the outcome of this procedure before publication. 
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Figure 20: Double-differential distribution with respect to the dijet invariant mass M 12 
and y* at the Tevatron. In the absolute prediction (uppermost plot) the cross section is 
divided by the bin width in y*. 
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Figure 21: Double-differential distribution with respect to the transverse momentum of 
the leading jet /c T1 and y* at the Tevatron. In the absolute prediction (uppermost plot) 
the cross section is divided by the bin width in y*. 
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Figure 22: Comparison of the weak one- loop correction (^ e °£ p to the transverse- momentum 
spectrum of the leading jet as obtained from our calculation to the result of Ref. 



For the integrated cross section with minimal cuts we find that the weak corrections are 
negligible, typically staying below the per-cent level, both at the LHC and the Tevatron. 
However, the electroweak Sudakov logarithms affect the tails of kinematic distributions 
that are sensitive to the high energy scales of the hard scattering process. We have 
discussed this feature in some detail, considering the distributions in the dijet invariant 
mass (M 12 ), and in the transverse momenta of the leading (/ct,i) and the subleading (&t,2) 
jets. 

For the \fs = 14 TeV LHC setup, we observe weak loop corrections of O {ot 2 a) reaching 
up to —12% (—16%) for a transverse momentum of k^ = 3 TeV of the leading (subleading) 
jet, whereas the dijet invariant-mass distribution only receives weak corrections up to —6% 
for M12 = 6 TeV. This difference is explained by the fact that observables based on specific 
ranges in M\2 are not dominated by the Sudakov regime (large energies at fixed angles) 
for large M 12 values, but rather characterized by the Regge (forward) regime. The weak 
corrections to dijet production at the Tevatron show similar features, though their size 
is smaller (—10% for k^ = 800 GeV) than for the LHC because of the smaller scattering 
energy. 

The LO EW contributions of O (a s a, a 2 ) turn out to be of the same order of mag- 
nitude as the weak loop corrections. At the LHC, these two types of corrections partially 
cancel, but the degree of this cancellation depends on the chosen observable, setup, and 
cuts, so that the full calculation is necessary in order to correctly include the considered 
electroweak effects. At the Tevatron the LO electroweak corrections are somewhat smaller 
than the weak loop effects, but of the same sign at high transverse jet momenta, thus 
somewhat enhancing the electroweak effects to —12% at k^ = 800 GeV. 

The electroweak corrections considered in this paper are supposed to be the by far 
dominant electroweak effects in dijet production at hadron colliders. Their numerical 
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impact of 10—20% in the TeV range is not negligible and will certainly play a significant 
role once the NNLO QCD corrections are known. In contrast to the weak corrections, 
which get dominated by the large Sudakov logarithms at high energies, the so-far-neglected 
photonic loop corrections do not receive particular enhancements over their parametric 
suppression by the electromagnetic coupling a. The calculation of these effects, which are 
expected to stay at the few-per-cent level, are left to the future. 
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